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PROBLEMS FOR ARBITRARY CASGADE3. OF AIRFOILS 
By William Mutterperl 



SUMMARY 

Solutions of the direct potential problem for a 
cascade - namely, to find the potential flow past an 
arbitrary cascade of airfoils - and of the Inverse prob- 
lem - namely, to find an airfoil having a prescribed 
velocity distribution In cascade - are presented. The 
methods used represent the extension to cascades of the 
Cartesian mapping function method of conformal trans- 
formation. Numerical examples of the direct and Inverse 
methods are given. The numerical labor required Is con- 
servatively estimated as double that Involved for iso- 
lated airfoils by the corresponding methods. 

INTRODUCTION 

The potential flow throu(7h a cascade of airfoils 
bears approximately the sarre relation to the real flow 
through axial and through some types of centrifugal 
turbines and compressors that the potential flow about 
Isolated airfoil sections bears to the real flow about 
airplane w5.ngs. A knowledge of the potential flow 
through a cascade is therefore recognized as basically 
Important for the aerodynamic study and design of such 
machines. The relative difficulty of making measure- 
ments on airfoils In cascade further accentuates the 
necessity of a potential-flow solution. 

As was the case for Isolated airfoils, there are two 
potential problems for the cascade that are of primary 
practical Importance. The first, or direct problem. Is 
that of determining the potential flow past a given arbi- 
trary cascade of airfoils. The second, or Inverse prob- 
lem. Is that of deriving an alrfoll section to have a 
prescribed surface pressure distribution In cascade. 

The theoretical methods available for the solution 
of these problems are, roughly speaking, of three kinds, 
namely: 
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(1) Methods that regapd a blade of the cascade as 
an isolated airfoil operating In a flow composed of a 
free-stream velocity and a disturbance velocity due to 
all the other blades of the cascade. The method applies 
best to cascades of thin airfoils with small solidity. 
(See reference 1, p. 70, and bibliography contained- 
therein.) 

(2) Stream-filament methods, which regard the 
space between the blades of the cascade as channels of 
varying area but in which the streamlines are uniform 
or of simple curvature. These methods apply best to 
cascades of high solidity in which, moreover, the flow 
is smooth (shock free) at entrance. 

(5) Methods based on conformal transformation of the 
cascade. These methods may be subdivided as follows; 

(a) Methods based on the concept of the equiva- 
lent cascade of flat plates; that is, the cascade 

of flat plates with spacing equal to that of the 
given cascade, with blade angle equal to the zero- 
lift angle of the given cascade, and into which the 
given cascade can be transformed conformally. 
Rjctenalve use of this concept la made in reference 1, 
on the basis of which are fiven appro:clp.ate solu- 
tions of the direct and inverse problems for cascades 
of various types of shape and for various ranges of 
solidity. The solutions are approximate mainly 
because of the methods given for the determination 
of the equivalent cascade from the given cascade 
or vice versa. 

(b) Particular conformal transformations that 
yield special classes of shape for which the flow 
can be calculated exactly (such as those of refer- 
ence 1, p. £5, and reference 2) . 

(c) Sxact methods for arbitrary airfoils 
or pressure distributions in cascade, such as 
those that exist for isolated airfoils. For 
this purpose, Welnig (reference 1, p. 90) uti- 
lizes the basic known transformation from a cas- 
cade of flat plates to a single circle. By this 
transformation, the given cascade transforms to 

a near circle. The near circle, p» -plane, is then 
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transformed to a circle, p-plane, by a modification of 
the Theodoraen-Garrlck procediire . (xef.erfxice JS) Qonal sting 
In the use of log dp '/dp Instead of log p'/p as 
the mapping function. The transformation of the given 
cascade to the near circle appears to be much more 
laborious and the near circle much more different from 
a circle than is the case for the transformation of 
isolated airfoils. This general procedure has recently 
been treated in reference 4. A similar method has also 
been used in a British paper by A. R, Howell. The 
basic transformation that reduces the given cascade to 
a single shape is taken as t = tanh' z, where z is the 
physical plane. Inasmuch as the resulting f -plane shape 
is not of near-circle type, several Joukowski transforma- 
tions are applied to produce a near circle. The transfor- 
mation from near circle to circle is then accomplished by 
the Theodorsen-Garrick metho('. This method, too, involves 
an excessive amount of numerical labor. 

In the present paper, the Cartesian mapping function 
method of reference 5 is extended to the solution of the 
direct and inverse potential problems for arbitrary cas- 
cades of airfoils. After an exposition of igctne basic 
properties of the Cartesian mapping function for cascades, 
procedures are given for the solution of the direct and 
inverse problems. The procedures are then Illustrated 
by numerical example 9. 

Acknowledgn-ent Is made to Mr?, Lois Evans Loran the 
computing staff of the lan^^ley ivill~»^cale tunnel for her 
assistance in making the calculations. 

SYMBOLS 

z = X + iy plane of cascade of airfoils 

£ = ^ + iT] plane of cascade of straight lines 



p plane of imit circle 

q> central angle of circle 

^x component of Cartesian mapping function (CMF) 

parallel to chord 

Ay component of Cartesian mapping function per- 

pendicular to chord 

T displacement constant for locating airfoil 
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blade aiigle of cascade, measured between 
chord line and nomal to cascade direction 

parameter of transformation from cascade of 
straight lines to unit circle 

central angle of unit circle that corresponds 
to extremity of flat plate in cascade 
(equation (7)) 

chord 

solidity of cascade of straight lines 
\/co8^3 + sinh^K (equation (12)) ■ 

adjustnent constant (equations (13) to (15)) 
velocity at airfoil surface 
velocity at surface of unit circle 
relative mean velocity in physical plane 
complex potential in circle p-plane 
circulation 

relative stream angle (See fig. 5.) 
tan Xj^ " ^ tanh K (equation (27)) 

corresponds to ideal stream anfrle (equa- 
tion (32)) 

tan (pp 

tan (fvpj^ = ^g^^ (See equation (21) ) 

lift force per unit lencth on airfoil in 
cascade, perpendicular to relative mean 
velocity 

density 



lift coefficient 
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a angle of attack, meapured between mean rela- 

' *' ■ ■ ■' tlve stream dlre&tlon-and chord line 

pm zero lift angle, angle of ^attack for zero 

lift 

0 ■ angle of deviation of relative flow from 

leading- or trai ling-edge direction of 
airfoil (See fig. 5.0 

k denominator of equation (30.) 

P pressure coefficient jl - ([^}^J 

Subscripts; 

N leading edge (nose) 

T trailing edge 

u upstream 
d downstream 

APPLICATION OP CARTESIAN TRAPPING FUNCTION TO CASCADES 



The general scheme of application of the ClIF to 
cascades of airfoils was indicated in reference 5. 
Basically, the method consists in relating conformally 
the points of an airfoil in cascade to the points of 
its extended chord line by means of the vector differ- 
ence between pairs of corresponding points . This 
vector difference, called the Cartesian mapping function 
(CMF), is a function of position in any of the mapping 
planes used. It is regular and periodic in the entire 
region outside the airfoil boundaries, in the corre- 
sponding region outside the straight lines, and is 
regular in the entire region outside the single unit 
circle into which the preceding regions can be confor- 
mally transformed. Instead of the simple Joukowskl 
transformation from the Isolated straight line to the 
circle, however, as was used for* the Isolated airfoil, 
the more complicated but well-known transformation from 



6 



NACA ARR No. L4K22b 



a cascade of straight lines to a clz*cle Is now necessary 
to express the CUF as a regular function outside the 
circle. Various preparatory details of these tranaforma- 
tlons are next set down, prior to a discussion of the 
general direct and Inverse problems for cascades. 

CMF for cascades .- The CMP from the cascade of alr- 
folls to the cascade of flat plates Is taken as z - £. 
(See fig. 1.) Because z - i Is regular everywhere 
outside the cascade of flat plates - In particular. In 
the Infinite strip outside one plate that maps Into the 
entire region outside the \anlt circle - It Is regular 
everywhere outside the circle and is therefore expres- 
sible by the Inverse power series^ 

On the corresponding boundaries of airfoil, flat plate, 
and circle 

z - £ = (x - §) + ly 

S Ax + lAy (2a) 

p = el<P (2b) 



and 



Substitution of equations (2a) and (2b) in equation (1), 
with Cn = aji + Ibni gives for Ax and Ay the con- 
Jugate Fourier series 

Ax(q)) = + ^^^n "•" ^^^n sin np (3) 



■'■The mapping function chosen is periodic with respect 
to the same strip in the £ -plane that maps into the 
entire p-plane region outside the unit circle. The 
necessity of introducing a cut in the infinitely many 
sheeted p-plane is thereby eliminated. 
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ea ■ ... ■ cw 

Ay(<p) = bo + ^]^n " ^ ftp sin n(p 



(4) 



The corresponding Integral relations between Ax and 
Ay are the sanie as In reference 5« 

Flat-plate transformation .- The transformation 

from the cascade of flat plates, S-plane, to a single 
unit circle, p-plane, which makes the CMP a function of 
the central angle qp of the circle. Is (fig. 2 and 
references 1 and 6) 

S - T = e-lP loge 4^ + l°fie <5) 

e"- - p p - e""- 



This transformation Is periodic In J with a cascade 
spacing (period) of 2rr, The points t = +» corre- 
spond to p = ±e^. The constant displacement t Is 
inserted for future convenience In locating the cascade 
of airfoils. The correspondence between points on the 
flat plate and points on the circle is obtained by 

substituting L = S + P = equation (5). 

The result Is 



£ = T + cos p logg cosh K + COS CP ^3 ^^^^^-isin^ 

® cosh K- cos <p slnh K 

71= 0 



(6) 



The suiples Pj^, ^ that correspond to the extreml* 

ties of the flat plate are obtained by maximizing g In 
equation (6) with respect to 9* The resulting condi- 
tion Is 

tan Pj^ = tan p tanh E (7) 

Substitution of equation (7) in eq\aation (6) gives the 
.„«aity = &SSig^ of a coaae of .tralght 
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lines In tema of the Islade angle p and parameter K 
as 



COS p loge + cos^p + cos p 



slnh K 



+ sin p tan"^ P \ (8) 

V slnh^K + cos^p J 



It Is remarked that. In the applications which are 
to be given, values of £ corresponding to selected 
cp-values will be determined from equation (6). The ■ 
relatively more difficult Inverse calculations of qj 
for a given or the still more difficult calcula- 

tion of p for a given i by equation (5), will not 
be necessary. Therein lies the essential numerical 
simplicity of the CW method for cascades. 

Airfoil position and adjustments In terns of CMF .- 
The coordinates x, y of the airfoil corresponding to 
the points at the angular positions flp on the unit 
circle are obtained from equations (2a) and (6) , as 



X = T + cos p loge K + cos <P 

cosh K - cos q) 



+ 2 sin p tan-1 -—tT^ + Ax(<p) (9) 

Sinn K 



y = Ay((p) (10) 



As was the case for Isolated airfoils, it will be 
necessary to find for cascades the points on the circle 
corresponding to the chordwlse extremities of the air- 
foil. The abscissa x In equation (9) is therefore 
maximized with respect to 9. The resulting condition 
is 
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6Ax _ C (cp - ^^^^ 



slnh^K + 'aln^qp 
where 



0=2 y ooB^p + slnh^K (12) 



and la defined by equation (7). Equation (11) la 

solved (graphically) for the angles <Pjj and qprp corre^ 

spondlng to the leading and trailing edges of the air- 
foil (defined as the extremities of the airfoil abscis- 
sas) . 

The solution of equation (11) will usually be 
carried out In connection with such an adjustment of the 
airfoil that the extremities are symmetrically located 
and add up to a desired chord (solidity). The adjust- 
ments are similar to the horlzsontal and vertical 
stretchings for the Isolated airfoil. The horizontal 
stretching adjustment consists In finding the length of 
straight line, that la, K, by equation (8) since |3 
Is the known angle of the airfoil chord line, which with 
a given Ax(q>) will place the airfoil extremities at 

+^ = ±TTa, where o la the given solidity and c is 

the airfoil chord. The vertical stretching adjustment 
consists In finding, for a given length and angle of 
straight line, the factor S by which to multiply a 
known Ax (9) In order that the resulting airfoil will 
have Its extremities at ^ = ±Tr<7, where again a Is 

the desired solidity. The equations for these adjust- 
ments are obtained by equating the abscissas of the 
airfoil extremities to their desired values; that Is, 

^m - ^- ^° 

cosh K + cos cpLj 

= T + cos p log_ r-rr . 

cosh K - .cos (jj^j 

+2 Bin p tan-1 + SAXj, (13) 
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Xj, = -| = -TTO 

cosh E + cos 9q, 

= T + COS p logfl r~= 

°e cosh K - COS qvp 

^ sin Q)m 

+-2 Sin fl tan"-"- , ^ 1 + SAxr, (14) 

slnh K 



In addition, equation (11) for the airfoil extremities 
Is rewritten as 

„dAx ^ C am (cp - ^^gj 
Blnh^K + sln^qp 



Subtraction of equation (14) from equation (13) 
yields 

(cosh K + cos (PjjN /cosh K - cos aw,\ 
i- ) ( ^ ) 
cosh K - cos qpjj / Vcosh K + cos flpTy 

.-2 sin pjta.-l(I^)- tan-l(l^)] 

+ S(Ax^ - Axj) (16) 

Addition of equations (13) and (14) yields 

(cosh K + cos flPiA / cosh K + cos (OA 
) ( ) 
cosh K - cos cpv/ \ cosh K - cos avn/ 

+ 2 sin P tan-l(-^^;^?r~)+ tan'^ r]"" 

|_ \alnh K/ \alnh K/J 



+-s(Axjj + AXfp) 



(17) 
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In the horlzont&l stretching adjustment, equations (15) 
and (16) with 3=1 ard solved simultaneously for 
cpjjV and --(P|j,>' - wl4:b'- 3 . and . . a. ., iSlven. In the vertical 

stretching adjustment, equations (15) and (16) are solved 
simultaneously for S, o^, and (Prp, vlth E, 0, and 

a given. Finally, equation (17) Is solved for t. 

Velocity and pressure distribution's .- Once the oon- 
formal correspondence between a cascade of airfoils and 
the unit circle has been obtained - that Is Ax((P) , 
^y(flp) » K, op^, and cp^j - the velocity at the airfoil 

surface Is given as usual by the product of the velocity 
at the corresponding point of the circle and the 
stretching factor from the circle to the airfoil, as 



ap 

The velocity In the circle plane. le that due to a 
superposition of sources and vortlce? at the points 

p = ±e (fig. 3) , which reproduce the de«!ired flow 
conditions at ±oo In the physical plane and which 
maintain the unit circle, and hence the cascade of air- 
foils, as streamlines (reference 6) . The effect of a 

source at p = e^ Is to produce a velocity at z - a 

frftm rlfTht to left of rfia^nitude that 1p, the 

co"iplex velocity Is -^c^*'* where m Is the total flux 

of the source. This fact is evident from the complex 
potential for a source and the transformation (equa- 
tion (5)) between p-plane and ^-plane. The additional 
effect of an equal source at the image point p = e"^, 
together with an equal sink at the origin. Is to make 
the circle (hence the cascade of flat plates and air- 
foil?) a streamline while maintaining the same flow 
condition at Infinity. Similarly, the effect of a 
counterclockwise vortex at p = e^ Is to produce a 

downward velocity at a =- « of magnitude that Is, 

the complex velocity is ^e^P, where k Is the vortex 
circulation. The additional effect of an equal and 

rr 

opposite vortex at p = e ^ is to make the circle a 
streamline. 



If the flow conditions at infinity in the physical 
plane are taken as consisting of a mean velocity V at 
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GUigle \ (fig* 1), on which Is superimposed the tan- 
gential velocity at z = to correspond to the 

circulation about the airfoil r (positive clockwise), 
the complex potential in the circle plane is that due to 
the system of sources, sinks, and vortices indicated in 
figure 3; namely, 

W = -V cos X log|i^^)^-t^^ 

-IV sin \ logfe-^I- ^P + 

Vp + e^/\p - e~V 



The circulation T is fixed by the Kutta condition of 
smooth flov at the trailing edge in the physical plane 

or dV//dp = 0 at p = e^^ In equation (19) . The 
resulting circulation is 



Jsin^qjT, + sinh^K 

r = 4TrV - ^ ^ ; : sin C^K - ^) (20) 

slnh K cosh K / 



where 



tan flPbi 



It is noted here that the resultant force on an 
airfoil of the cascade is, by elementary cascade theory. 



I = pvr 



(22) 
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where I acts perpendicular to the relative loeah 
velocity .. V. ...The. lift coefficient of the section based 
on the relative meem velocity V Is" therefore, with 
c = 2va, 



Setting 



and 



\ B a + p (24a) 



qVpK = IT - + P (24b) 



and from equations (20) and (23), the lift coefficient 
of an arbitrary airfoil In cascade Is 



4 \|Bln2<p^ + alnh^K 
°l = 5 slnh K cosh K + (25) 



The angles a and Pip are, respectively, the mean 

relative stream angle of attack and the zero-lift angle 
measured from the airfoil chord line. Equation (25) 
Is seen to be similar in form to the corresponding 
equation for isolated airfoils, to which It reduces In 
the limit of zero solidity. 

With the circulation of equation (20), the velocity 
distribution on the unit circle is obtained by differ- 
entiation of equation (19) and substitution of p = e^*"'. 
The result Is 



Vp(CP) 



= 2i 



oos^X t -^^"^ [sm ((P-Xj5).sln (cp^ - ^k}] 
sln^qp + slnb^K ^ n ^ xv/j 



(26) 
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where 

tan = tan \ tazih K (27) 

The stretching factor- j-^j from the circle plane 

to the flat-plate cascade plane Is obtained by differ- 
entiation of equation (5) . The result Is equal to 
Vp((p) (equation (26)) with the circulation term omitted 

and \ = p. (The transformation (5) can. In fact, be 
derived (reference 6) by a hydrodynamlcal argument 
similar to that leading to equation (19).) Thus, 



p" '" o Sin (q) - pt^ 
cos^p + slnh^K ^ (28) 

sln'^cp + slnh'^K 



where has been defined by equation (7) , and 



depending on the sign of cp - pj^. 



The over-all stretching factor j-^ j from the 

circle plane to the physical plane Is obtained by dif- 
ferentiation of equation (2a), 



z - £ = Ax + lAy 

_ d£ . 

- r^li«-i<P 4. ■ . ■ dAy 



op dp dp 



or 



NACA ARR No. L4K22b 



15 



The velocity distribution at the surface of an 
arbitrary airfoil In cascade Is therefore, by equa- 
tions (18) an* =(S6.),J:o.,( 2a) » ■ 



2 \|oos^\ + alnh^K[8ln (<p - Xj;) - sin (cpkp - Xj^jl 
(sm^. . slnh%) . . 



T ' 



(SO) 



The nianerator of this expression Is seen to be dependent 
on, but the denominator Independent of, the stream 
angle \. 

The pressure coefficient Is given by Bernoulli's 
equation 



P = 1 - K^J (31) 



where P Is the ratio of the local increment of static 
pressure over the pressure at which the velocity Is the 

mean stream value V to the dynamic pressure ijPV^. 

Integration of this pressure coefficient aroxind the air- 
foil section yields the lift coefficient defined by 
equation (23) . 

The expression for the ideal angle of attack may 
be noted. This angle is the mean relative stream angle 
for which the leading edge of the airfoil, corresponding 
to (p = % on the unit circle, is a stagnation point. 

It is obtained from the condition Vp((p) =0 In equa- 
tion (26) as 



(The corresponding angle of attack measured from the 
"x-axis" chord line is given by equation (24a),) 
Conversely, for any stream angle X, the front stagna- 
tion point is given by 



cpjj = w + 2X^ - qpip 



(33) 
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Superposition of solutlona .- Equations (3) and (4) 
show that the sum of two CMP 'a Is itself a CMP. Thus, 
as for isolated airfoils, the superposition of solutions 
for airfoils in cascade Is permlsslbie. Various types 
of superposition are possible because the solidity and 
stagger parameters E and p used with a resultant 
CMP need not be the same as either of these parameters 
used with the component CMP*s» The use of superposition 
is Implicit in the horizontal and vertical adjustment 
derived previously for an airfoil in cascade. 

Equivalent flat-plate cascade The many useful 
results of reference 1, which are derived and presented 
on the conceptual basis of the equivalent cascade of 
flat plates, make it worth while to note some of the 
properties of the equivalent cascade. In general, the 
conformal transformation of a cascade of airfoils into 
a unit circle determines at the same time every possible 
cascade of flat plates into which the given cascade can 
be confomally transformed, with the spaces at +ea cor- 
responding and equal. All such cascades of flat plates 
transform to the unit circle with the same locations for 

the points p = ±eK corresponding to the points 
in the cascade planes. The parameter K of the flat- 
plate transformation (equation (5)) is therefore fixed. 
Of the infinity of flat-plate cascades corresponding 
to the remaining parameter, namely, the stagger angle p, 
that one of which the stagger angle equals the zero-lift 
stream angle of the given cascade is by definition the 
equivalent cascade of flat plates. This stagger angle 
is (D^j^ by equation (20) . Equation (21) then shows 

that the angle (prp corresponding to the trailing edge 

of the airfoil in cascade corresponds also to the 
trailing edge of the equivalent flat plate. 

The solidity of the equivalent flat-plate cascade, 
as given by equation (8), is evidently somewhat larger 
than that of the given cascade. For a given stream 
angle \, equation (20) shows that the circulation T 
and thus, by equation (23), the product ac^ are the 

same for the two cascades. 

The velocity distribution in the airfoil cascade 
at one stream angle may be converted to a velocity 
distribution at any other stream angle by a knowledge 
of the equivalent cascade, because the ratio of the 
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velocities at two different stream angles Is the ratio 
of . the. cprres ponding ntunerators In equation (30). 
Inasmuch as the" constants "of ^the -numerator In equa-.. . 
tlon (30) are the same for the- given cascade and the 
equivalent cascade « the desired velocity ratio equals 
that for the equivalent cascade. This result Is 
analogous to that for Isolated airfoils and their 
equivalent flat plates. It Is useful when the equiva- 
lent cascade can be easily determined to a satisfactory 
approximation (cf . reference 1) . If the corre-spondence 
between an airfoil In cascade and Its equivalent flat 
plate In cascade Is known with sufficient accuracy, 
however, the CMP of the transformation Ax((p) + lAy(<p) 
can be set up and the velocity distribution obtained 
directly by equation (30) . 



SOLUTION OF DIRECT POTENTIAL PROBLEL' FOR CASCADES 



The direct potential problem for cascades of air- 
foils Is that of determining the general potential flow 
past, and in particular the velocity at the surface of, 
a known arbitrary airfoil In a cascade of given 
solidity a and stagger angle p (fig- The solu- 

tion Is analogous to that In reference 5 for Isolated 
airfoils. The method Is one of si^ccesslve approxima- 
tion, whereby an airfoil of which the CI.:P transformation 
to Its "chord-line" cascade and hence to the unit circle 
Is known Is compared --vlth the given airfoil to obtain 
an IncreF-ental CMP. The incremental Clv'P plus the known 
CMP yields an airfoil that may again be compared with 
the given airfoil, and so on. 

The steps in the process are as follows: 

Step 1: An airfoil in a cascade of stagger p and 

solidity a la given and an initial airfoil is assumed, 

for which are known the CMP AXq(cp), Ayf,(cp), the 

angles q)w , <pm corresponding to the leading and 
^•0 -^o 

trailing edges, the cascade parameter K^, and the 
translation constant - all adjusted for a cascade of 

the given stagger and solidity. The chordwlse locations 
Xq for this airfoil corresponding to a set of evenly 

spaced values of (p on the unit circle are calculated 
by equation (9) . 
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Step 2; At tbs ohordwlse locations Xq, the dif- 
ferences 6jj = Ay^ - Ay© between the ordinate 3 of the 

given airfoil and those of the assumed airfoil are 
obtained, graphically or otherwise. It is Important that 
the given airfoil be located relative to the coordinate 
axes as shown in figure 1 and that, when the airfoil 
ordlnates Ay^ are obtained, the airfoil boundary be 

traversed (counterclockwise) in the sare sense as Is the 
extended chord line . 

Step 3s The corresponding conjugate chordwise- 
dlsplacement lncrementf< 6x]l calculated In 

accordance with equations (3) and (4) . The slopes 
d5x]^/dcp and dSy^/dcp are also computed by means of 

the corresponding derivative Fourier series. Details 
of this calculation are given In appendix A of refer- 
ence 5. 

Step 4s The f irf?t-approxiraation CMP, namely, 
Axj = AXq + ^^1 " ^^o ^yi» corre- 

sponding slopes dAx-j^/d<p, dAy^^/dcp are now known. In 
conjunction with the constants Kq, (Pjt , qptp , and 

of the initial airfoil, this C'fF would, in general, 
determine an airfoil (equations (9) and (10)) of which 
the extremities (equations (11), (13), and (14)) would 
not add up to the desired chord (c = Stto) , although the 
stagger angle p (of the x-axls chord) is correct. 
The constants Kq, 'PNq* *PTq# ^^'^ """o therefore 

corrected by a horizontal stretching adjustment (equa- 
tions (15), (16), and (17)) such that the resulting 
constants K^j %i » ^-tt "^1 yield an airfoil of 

the desired chord. 

Step 5: The new chordwlse locations x-j^ are cal- 
culated by equation (9) with the constant? derived in 
step 4. The resulting airfoil, of coordinates x^ , Ay^^, 
is compered with the given airfoil. Tf the agreement 
is not satisfactorily close, steps 2 to 5 are repeated, 
all subscripts being advanced by one. 

Step 6: If the agreement is close enough, the exact 
velocity distribution of the airfoil just derived (which 
is thus a close enough approximation to the velocity dis- 
tribution of the given airfoil) is calculated as indicated 
in the section "Velocity and Pressure Distributions." 
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SOLUTION OP INVERSE POTENTIAL PROBLEM FOR CASCADES 



As compared with the Isolated airfoil, there are 
several possible variations of the Inverse problem for 
cascadesj depending on the over-all cascade data given 
In addition to the prescribed velocity distribution. 
These cascade data Include the ultlicate upstream and 
downstream flow angles considered as one quantity, the 
solidity a, and the stagger angle p. Of these three 
quantities. It appears that two may be specified 
Independently along with the prescribed velocity dis- 
tribution. 

The Inverse problem for cascades can be solved by 
the method of comparison of surface potentials (riven in 
reference 5 for the case of the isolated airfoil. A 
numerical example, howevfjr, disclosed several difficul- 
ties of application. The over-all adjustment for one- 
to-one correspondence of the true and approximate 
potential curve? was undu]y laborloixs. Furthermore, 
the first exact velocity distribution obtained oscil- 
lated about thp prescribed distribution through wide 
limits. This oscillation appeared to be caused by the 
somewhat indirect method of determining 6x((p), 
resultinp in inexact and oscillatory value." for the 
derivatives dAx/dcp, dAy/dcp, 

The more Ptrairhtf onivard method of derivatives, 
namely, determining dAx/dcp, dAy/dCD directly from the 
prescribed velocity distribution, was therefore resorted 
to. As remarlced in reference 5, the method can be 

regarded as based on the use of the function Ip ^^^ " 



This function is single-valued and regular everywhere 
outside the unit circle; It approaches the limit zero 
as p-^oo and on the unit circle Itself p = e^-^P 



tlves can be so determined as to satisfy the expres- 
sion for the prescribed velocity distribution (equa- 
tion (30) ) , the airfoil Itself can then be obtained 
by a slznple Integration. 

The determination of the derivatives can be made 
by successive approximation. The constants \, F are 
fixed to satisfy approximately the over-all cascade 



dp 



reduces to 




If, therefore, these derlva- 
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data previously mentioned. Inasmuch as the derivatives 
dAx/dop, dAy/dqp are of the same order of magnitude for 
the shapes to be dealt with. It is apparent from the 
structure of the radical in equation (30) that d<ix/d(p 
will exert a greater influence on the velocities, over 
moat of the airfoil than will dAy/dcp. To a given 
approximation, dfiix/d(p can therefore he solved for in 
terms of the prescribed velocity Vg and the data 

dAy/dcp, \, and K of the previous approximation. A 
new dAy/dcp can then be calculated as the function 
conjugate to this dAx/dcpj and \, K can be corrected 
to serve ap the basis for the next approximation. 

The steps of the inverse process are outlined as 
follows; 

Step 1: A velocity distribution expressed as a 
function of percentage chord or arc and the appropriate 
cascade data are given. The definition of the chord and 
the prescribed velocity distribution are taken as in 
reference 5. An initial airfoil is assumed, which ic 
Judged to produce approximately the desired conditions 
and for which the cascade CMP Ax^, AyQ is known. 

3tep 2: If the given cascade data are the solidity a 
and stagfer angle p, a horizontal adjustment la applied 
to the CLIP Axq, Aye to achieve these values; that is, 
equations (15) and (16) with S = 1 are solved for K^, 
cpjj^, and C!>r^ with the given values of a, 3. Thence, 
equation (17) is solved for Tq. 

Step 3: The lift coefficient cjq corresponding 

to the prescribed velocity distribution la obtained 
approximately a? the area under the curve of chordwlse 
pressure distribution. A correction to this c^, if 

necessary, is obtained by integrating the chordwlse 
perpendicular pressure components of the prescribed 
distribution on the initial airfoil.^ The required 



'Equations (22) and (23) are dynamical equations, holding 
only for simply connected airfoils to which Bernoulli's 
equation applies. In order to preserve the consistency 
of the CMP equations for figure-eight contours, the 
circulation r can be approximated and equation (20) 
can be used directly. 



NACA AHR No. L4E22b 



81 



atream angle \q for this cj^ and o la calculated 
from equations (23)- and- (20) by use of qpip^ and, . Kq. ^..pf 

step 2. If the ultimate upstream and downstream flow 
angles are specified Instead of the solidity^ an e::ctra 
relation between T and X Is known. If this relation 
and equation (23) are substituted In eqviatlon (20)^ an 
equation relating (p^, K, emd a results. Thld equa- 
tion la solved simultaneously with equations (15) and (16) 
for' <PNq* *'^o' These values then deter- 

mine Tq, Xqj and o^, A similar situation exists If 

the upstream and downstream angles and the solidity a. 
Instead of the stagger angle ' are specified. 

Step 4: 'With the constants derived, the chordwlse 
locations Xq are calculated by equation (9) for a set 
of evenly spaced cp-values. 

Step 5: The first approximate derivative dAx^/dq) 

Is calculated by solution of equation (30) with the 
constants already derived, the prescribed velocity Vg/V 

corresponding to the chordwlse locations x^, and the 

Initial dAyo/d(p. 

Thus, 

^ = [i] * i^J - (^T 

Two possible seta of values for dAx^^/dqp result from 

this calculation, depending on the sign of the square 
root In equation (34). It appears that for about one 
half the (p-values, one set of dAx^/dq) roots should be 

taken and for the other half, the other set. The 
resulting dAx^^/dop should be of the right order of 
magnitude to lead to a real airfoil. Any Imaginary 
roots that occur are replaced by the most favorable 
real values . (The assumption Is that a real, solution, 
exists and that any Imaginary values of dAx/dcp In 
equation (34) are the result of too poor an Initial 
choice of dAy/dqp-.) 
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Step 6: Ttie derivative dAy^/dcp conjugate to the 
dAX]^/d(p of step 5 le calcvilated In accordance with the 
derivatives of equations (3) and (4). The functions 
Ax;Lt*P)» Ay^Ccp) are ther. obtained by Integration. If 
the difference 

dSxi dAxj dA^Q 
dcp d(p dop 

Is small, It may be more convenient, because of smaller 
required accuracy, to compute the Increments dSy-j^/dcp, 

fiXj, and By^ and add them to the corresponding values 
of the previous approximation. 

Step 7: The horizontal adjustment of step 2 and 
the calculation of the stream anfLle X-^* step 3, Is made 

by using the f Irst-approxlnatlon CMF derived In step 6. 

Step 8 1 The chordwlse locations X]^ are calcu- 
lated as In step 4 by using the constants derived In 
step 7. 

Step 9: The first approximate airfoil and Its 
exact velocity distribution are calculated by equa- 
tions (9), (10), and (30). If the velocity distribution 
Is not satisfactorily close to that prescribed, the 
procedure Is continued with step 5, all subscripts being 
advanced by one. 



ILLUSTRATIVE EXAMPLES 



Direct method .- As an Illustration of the direct 
method, a cascade of solidity a = 1 and blade 
angle p = 45° was assumed, together with the given 
airfoil section shown In figure 4, of which the ordlnates 
are given In table I. In the resulting cascade arrange- 
ment (fig. 1) the flow Is from right to' left (compres- 
sion action) . 

As a preliminary step, the direct CMP method 
(reference 5) was applied to the given section considered 
acting as an Isolated airfoil. The resulting CMP and 
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transformation constants are given In tables IT and III, 
respectively. It Is remarked that the relatively poor 
Inl 1 1 al aipp r ox luist 1 on cho s en"j ■ nam. e ly / ' a J oiakow a kl al rf o 1 1 
of 6-percent camber and 12-percert thickness, necessi- 
tated three approximations before coincidence with the 
given airfoil (to a scale of chord length equal to 20 In.) 
was attained. 

For the airfoil In cascade, the procedure outlined 
In the section "Solution of Direct Potential Problem 
for Cascades" was followed. In accordance with step 1 
of the procedure, the Initial approximation was taken 
as the cascade of chord lines of unit solidity and 45 
blade annrle. The Initial CW was thus Axq = Ayg = 0. 

The parameter Kq, as determined from equation (8) 

with = 1 and 3 = 45°, was 0.3.'!0. The chordwlse 

locations Xq were calculated from equation (9), with 
T = Ax = 0. 

Steps 2 to 5 were then carried out. The resulting 
first approximate "airfoil" is shown in fif:ure 4. The 
abscissas are the chordwlse locations x-^Ccp) of the 
first approximation and the ordlnatec are those of the 
■flven airfoil at the chordwlse locations Xot'') ^^j® 
initial approximation, which are indicated 3n fip-ure 4. 
The first approximation is seen to have accomplished 
most of the required change from, the chord line to the 
given airfoil. The resulting ahape is not a physically 
real airfoil, however, because of the loop toward the 
trailing edge. 

A second approximation gave results ^hown in fig- 
ure 4. The agreement over most of the airfoil 1p good, 
although appreciable departure? over the front upper 
surface and a slight loop near the trailing edge still 
remained. A third approximation removed practically all 
the remaining discrepancies. The points that were -?tlll 
perceptibly different are indicated in flfure 4. A 
fourth approximation gave coincidence with the given air- 
foil to a scale of 20 Inches for the length of chord. 

The coincidence of the fourth approximate airfoil 
with the given airfoil was obtained for 24 points 
(besides the leading and trailing edges) . These points 
correspond to 24 evenly spaced values of (p on the unit 
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circle. It l3 seen from figure 4 that the airfoil 
points are sparsely distributed on the front upper sur- 
face and rear lower surface, as compared vrlth the rest 
of the airfoil. This distribution Is due to the 
Influence of the flat-plate transformation (equation (6)) 
In affectlnp the airfoil abscissas (equation (9)). 
Since the CMP was obtained In the form of 24- term Fourier 
expansions. It was thought that airfoil points obtained 
by Interpolation from the derived CMF might show some 
wavlness relative to the given airfoil In these regions. 
Five such Interpolated points are given In table II and 
are shown In figure 4. The wavlness appears to be 
negligible In this case. It Is noted that any desired 
spacing of points on the airfoil could be obtained by 
working with appropriate unevenly spaced cp-polnts on 
the unit circle; however, the determination of Ax((p) 
from Ay(cp) becomes more laborious. 

The slight wavlness In the Fourier expansions for 
Ax((p) and ^y(a)) gives rise to greater oscillations 
In the corresponding derivative expansions for dAx/dcp 
and dAy/dcp. In this example, the slopes were computed 
by harmonic synthesis of their Fourier expansions, and 
the results were smoothed out graphically to correspond 
to a faired curve through the 24 known values of ^x 
or Ay. The limits of oscillation of the derivative 
curves were small enouch that no appreciable error was 
believed to be Incurred In this process.^ 



In extreme cases, where the higher harmonics In the 
Ax{q)), Ay(cp) curves are large, the oscillations of 
dAx/d(p, dAy/dcp may lead to appreciable Inaccuracy of 
the velocity distribution. It Is possible, however, 
to solve Initially for the derivatives dAx/d(p, dAy/dip 
Instead of Ax(cp), Ay(cp) by working from the slopes 
dy/dx of the given airfoil, for 



d^ _ dy 

dx d(£ + Ax) 

_ dAy/dcp 
d^ dAx 
dqp dCD 
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The resulting CMP 'a and the transformation constants 
for the airfoil In cascade are given In tables II and III 
respectively^ --The - column of values of ■ k. unde^*- "Cascade* 
In table II gives the values of the denominator In equa- 
tion (30) for the velocity distribution. This factor 
does not change with. angle of Incidence of the flow. 

The variation of lift coefficient with angle of 
attack for the cascade arrangement (equation (25)) w^a 
obtained as 

ci = 4.65 sin (a.+ IOO4') 



where. It will be recalled, c^ Is based on the mean 

velocity V and pc Is the angle between the mean 
velocity V and the chosen chord line. The 11ft- 
coefflclent variation for the Isolated airfoil was 
obtained as 

ci = 6.86 sin (a + 11°22') 



. with the result that 



dAy _ dj/d^ dAx \ 
Bx\^dcp dcp / 

Thus, to a given approximation, dAy/dcp can be deter- 
mined from the known data of the previous approximation. 
To the given approximation, dAx/dflp is then the^ f\mc- 
tlon conjvigate to dAy/dcp. The derivative dAy/dcp can 
then be determined to a better approximation, etc. In 
each approximation, an exact airfoil is determined by 
Integration of dAx/d(p, dAy/dcp (the Integrated Fourier 
series are smoother than the original aeries) . This 
method for the solution of the direct problem is the 
counterpart of the method of derivatives in the Inverse 
problem. The accuracy of the derivatives dAx/dcp, dAy/d<p 
thus determined depends primarily on the accuracy to 
which the airfoil slopes dy/dx are known. Corre- 
sponding remarks apply to the isolated-airfoil case, In 
which, however, the higher harmonics are usually of 
less magnitude. 
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The ratio of lift-curve slopes la thus - g' - gg = 0.678 and 

the zero-lift angle for the cascade has heen reduced by 
1°18» . 

The equivalent cascade of flat plates has a stagger 
angle of 45° - 10^4' = 34°56» and the solidity, as 
determined by equation (8) with p = 34°56' and 
K = 0.2S8, Is Oi = 1.19. It may be noted that the 

ratio of lift-curve slope of the equivalent cascade to 
the Isolated-plate value Sir Is 0.624, vhlch Is about 
8 percent less than the value for the given cascade. 
The corresponding ratio for the chord-line cascade 
(p = 45°, a = 1) Is 0.844, or about 25 percent higher 
than that for the given cascade , If the equivalent 
cascade were based on the isolated-alrf oil data, namely, 
a, = 1.0912 and p = 45° - 11°22» = 33°38', the ratio 
of llft-curve slopes of the equivalent cascade Is 0.655, 
which Is a satisfactory approximation to the true value 
In this case. 



The variation with lift coefficient of the various 
flow angles of the cascade as well as the over-all 
preasure-rlge coefficient baaed on mean dynamic pressure 
Is given in table IV. The values were calculated by 
means of the velocity triangle of the cascade (fig. 5) . 
The deviations 5jf, 5ip of the upstream and downstream 

flow angles from the mean directions of the leading and 
trailing edges of the airfoil (fig. 5) are also given. 
The downstream flow angle \^ is seen to remain essen- 
tially constant as the upstream flow angle Is varied. 
In this respect, therefore, the cascade acts like one 
oi infinite solidity. The deviation Sip of the down- 
stream flow from the direction of the trailing edge Is 
large, however (about 18°) . 

The velocity distributions of both the Isolated 
airfoil and the airfoil in cascade, are shown for zero 
lift coefficient in figure 6 and for various other lift 
coefficients in figures 7 and 8. The velocity distri- 
butions at zero lift may be compared directly. Sxcept 
near the leading and trailing edges, the velocities for 
airfoils in cascade are higher than the velocities for 
the Isolated airfoil. This result may be ascribed to 
the constricting effect of two neighboring blades of 
the cascade on the streamlines between them. The 
velocity bump about 60 percent chord back on the upper 
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surface of the airfoil In cascade Is probably caused by 
the suddeji.rate of change of flow area In this region 
due to the jireserice bf the upper neighboring blade 
(fig. 5). At the leading edge, the velocity peak Is 
lower for the airfoil In cascade than for the Isolated 
airfoil, probably because of the effect of the lower 
adjacent airfoil In deflecting the upstream flow upward 
toward the nose. At the trailing edge, the velocity Is 
lower for the airfoil In cascade than for the Isolated 
airfoil, probably because the upper adjacent airfoil 
provides a divergent channel for the streamlines In 
this region. 

The proper basis ftor comparison of velocity distri- 
butions of the isolated and cascade airfoils at lift 
coefficients other than zero is not clear. Because 
the cascade velocities (flg« 8) have been expressed as 
fractions of the mean velocity, the velocities toward 
the leading edge are higher than if the ultimate up- 
stream velocity had been used as a base; whereas, tbwai'd ' 
the trailing edge, the velocities are lower than if the 
ultimate downstream velocities had been used as a base. 
The ultimate upstream and downstream velocities for the 
various lift coefficients are indicated by the hori- 
zontal lines in figure 8. In general. It appears that 
the adverse velocity gradients near the leading edge are 
less for the airfoil In cascade than for the Isolated 
airfoil, whereas the adverse velocity gradients near 
the upper-surface trailing edge are greater in the 
cascade case. 

The "ideal" velocity distributions, which produce 
a stagnation point at the chordwlse extremity of the 
airfoil, are those for isolated ~ ^'^^ ^^'^ 
^^cascade ~ It is seen that small velocity peaks 

are present on the lower-surface leading edge under 
these conditions. These peaks are removed in the 
isolated-airfoil case at c^ = 1.50 and for the cascade 

at about c^ = 1.30 without producing peaks on the 

upper surface. The incidence of the upstream flow rela- 
tive to the chord line is 15°24' in the case of 
Ol = 1.30. 

The four approximations in this exaiaple were carried 
out without any attempt at picking off the ordinates In 
order to anticipate the resulting change in chordwlse 
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location. It may be expected that, vhen more favorable 
Initial approximations than the chord line are avail- 
able, and with some experience In anticipating changes 
In chordwlse locations, no more than two approximations 
will be necessary for accurate result? - that Is, 
results less than 1 percent In error. 

Finally, It may be noted that the CMP derived for 
the given airfoil can be used as a good Initial approxi- 
mation not only for "Imllar airfoils \ti compressor 
action but also for airfoils In turbine action. For 
this latter application, the airfoil would be drawn 
with the camber on the same side of the chord line as 
In figure 1. The flow Is then from 3 eft to right; that 
Is, \ Is Increased by 180°. This procedure appears 
more favorable nur'.erlcally than drawing the airfoil with 
the camber on the other side of the chord line, which 
would maintain the flow from right to left. 

Tnverse method .- In order to illustrate the Inverse 
method for cascades, the velocity distribution at 
ci = 1.3 of the airfoil Just analyzed was modified as 

Indicated In figure 9. The corresponding prescribed 
pressure distribution is constant up to 50 percent 
chord and thereafter increases linearly to the trailing 
edge. The prescribed cascade data were taken as 
a = 1, p = 45 . The ultimate upstream, and downstream 
angles, or blade circulation, were expected to be about 
the same as those for the initial airfoil In cascade at 



The procedure already outlined was followed. The 
initial airfoil wan chosen as the one just analyzed by 
the direct m.ethod (fig. 10) . The C>.!F and associated 
constants of the initial airfoil are given in tables II 
and III. The Initial CT.?P wae actually taken as 1.145 
times the values of table II, because the use of the 
method of potentials had Indicated thi? factor as a 
preliminary over-all adjustment. The chordwlse loca- 
tions of table II were consequently corrected. The 
initial lift coefficient was taken as the area enclosed 



The first-approximation dAx^/dqp was thereupon 
calculated from equation (34) and plotted in figure 11. 
Except at (p = 15 , this curve represents the solution 



= 1.3. 



by the prescribed 




namely, ci 



= 1.26. 
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of equation (34) with the negative square root for 

jS.Q** it, (p % .180® and the solution with the positive root 

for the remaining q>-vatue8. At" cp i= 15°"; the calou— . . 

lated roots were complex. The -curve was therefore 

d<p 

faired through (p = 15 to complete the specification 
of dAx^/dtp. The increment 

dCx^ dAzi dAzg 
dtp . d(p dcp 

was then U5ied to calculate the first-approximation oon- 
Jufrate increment 

d6yi d&yi _ dAyp 
dcp ~ dcp " dtp 

and the Incremental airfoil CMP dx-^ + iSy]^. The 

Rungte 24-polnt schedule (f?ee appendix A of reference 5) 
was used for the^e calculations. The horizontal adjust- 
ment then applied to this CWP was such an to maintain 
0=1, Q = 45*^ (with reppect to the chord chosen for 
the Initial airfoil, the eo-called x-axl? chord) . The 
first apprDxlTiate CstP is listed In talile V and the 
associated constantn In table VT. The first approxi- 
mate airfoil is phown in figure 10 and It? exact 
velocity diptrlbutlon, for c^, = 1.26, in figure 9. 

The cascade flow anplea for the derived sections are 
given in table VII. 

The change? in velocity distribution and airfoil 
contour and position from their initial values are seen 
to be considerable. The mean line of the derived sec- 
tion Is reflexed toward the trailing edge, probably 
because of the influence on the "mean" streamline of the 
vorticlty toward- the leading edge on the upper neigh- 
boring airfoil (fig. 5). The maximum thickness of the 
derived section has been Increased from 10 percent to 
about 15 percent. Furthermore, the section as a whole 

1° 

has been rotated about 2-^ , with the result that the 

cascade blade angle, although equal to 45° as meas- 
ured from the initial or x-axis chord line, is about 
42,5° as measured from the "longest-line" chord in the 
derived airfoil. A slight change in cascade solidity 
is similarly brought about. 
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A second approximation to the prescribed velocity 
distribution was next carried out, with the Initially 
chosen chord maintained as the "chord" of the derived 
airfoils . It was evident that the poor first approxi- 
mation on the upper- surface leading edge, particularly 
for 9 = 15°, was due originally to fairing the 
d5x-i 

— g~-curve through this point too smoothly. The 

resulting value of dAx;j^/dcp was not the best possible 

real value, which la in this case the one that makes 
the velocity on the airfoil at this point as high as 
possible or, from equation (30), 



/dAx\ ^ [dn 



/dAxgX 

This value was therefore given to (— r—- I for the 

second approximation with the remaining values calcu- 
lated as in the first approximation. The incremental 
dGxo dAxp dAx-i 

derivative = —5— 5—= thus obtained la shown 

d(p dcp dcp 

in figure 11. The conjugate function dfiyg/dcp calcu- 
lated therefrom Is shovm in figure 11. The oscillatory 
nature of dOyg/dcp is caused by the large nagnltude of 

the higher harmonics in the Fourier series for dSxg/dcp. 

The Cyg-functl on, as obtained by synthesis of the 

Integrated Fourier aeries for d5yg/dcp, displayed very 

little oscillation. The slopes of the smoothly drawn 
6yg-curve were therefore measured graphically at all 

but the cp = 0°, 15° points. The ~^-curve thus 

obtained la seen from flfure 11 to be a riean curve 

dAyg 

through the calculated — ^^-valuep. This curve was 

used in the calculation of the second-approximation 
velocity distribution. Fse of the measured values of 
dSyg/dQ) as compared with the calculated values removed 

an oscillation of, at most, ±2 percent in the velocity 

distribution. Any error in measurement of the 
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-a lope a would therefore produce the quite negligible 
error of a email fraction of ±2 percent . 

The second approximation waa carried out for 
c^, = 1.26. The resulting CMF and associated constants 

are given in tables V and VI. The second-approximation 
airfoil and its exact velocity distribution are shoim 
in figures 10 and 9, respectively. The velocity distri- 
bution la seen to be converging satisfactorily to the 
prescribed values, though perhaps not so rapidly as 
would be the case for an isolated airfoil. The tendency 
toward formation of a velocity bump about 60 percent chord 
back on the upper surface is seen. As previously men- 
tioned for the direct method, this effect may be due to 
the sudden constriction of the stream caused by the 
forward part of the upper neighboring airfoil (fig. 5). 
The second approximate airfoil differs from the first 
approximate airfoil raalnly in having greater thickness 
toward the nose and leas over'-all camber. The greater 
thickness toward the nose is evidently the consequence 
of a higher required velocity on the upper surface at 
the loading edge. The reduced over-all comber appears 
to be caused by the necessity of a more constricted 
stream over the rear upper surface (higher required 
velocities) , The forward half of the upper neighboring ■ 
profile (fig. 5) shifts downward to accomplish this 
constriction. 



A sharp cusp, possibly even a slight, physically 
impossible loop, is evident at the trailing edge. 
Apparently, the prescribed velocity in this region la 
close to being unattainable with a physically real 
airfoil. The necessity of as sharp a cusp as possible 
at the trailing edge seems definite. 

Prom the practical point of view, two approxima- 
tions are considered sufficient in this case. The 
necessary modifications of the initial airfoil to 
produce the modified velocity distribution are clearly 
indicated. The results of this and the preceding 
example can now be used as good initial approximations 
in further modifications either of airfoil shape or of 
velocity distribution, or, to a lesser extent, of 
cascade geometry. 
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CONCLUSIONS 



1; The direct and Inverse Cartesian mapping 
fvmctlon (CMP) methods already used for Isolated air- 
foils can be applied directly and analogously to the 
cascade case. 

2. The cascade problem represents a more severe 
test of the CUP method than does the Isolated-alrfoll 
problem. Althouflth the convergence of the successive 
approximations Is slower In the case of cascades than 
In the corresponding Isolated-alrfoll cases, the 
numerical work for a. given accuracy Is considerably 
less than hitherto required, at least for cascade 
solidities near unity. 
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TABLE I.- AIRFOIL 0RDINATB3, DIRECT PROCESS 
[statlon8..a]3d^..oFdlnates 1a percent of airfoil chordj 



Upper surface 


Lower surface 


station 


Oralnate 


station 


Ordinate 


0 


0 


0 


0 


.23 


.97 


.77 


-.43 


.44 


1.24 


1.06 


-.44 


.88 


1.66 


1.62 


-.46 


2.03 


2.63 


2.97 


-.37 


4.43 


4.17 


5.57 


-.13 


6.87 


5.44 


8.13 


.18 


9.34 


6.56 


10.66 


.48 


14.34 


8.46 


15.66 


1.08 


19.38 


9.98 


20.62 


1.64 


24.44 


11.24 


25.56 


2.16 


29.55 


12.20 


30.45 


2.60 


34.64 


12.90 


35.36 


2.96 


39.72 


13.39 


40.28 


3.29 


44.82 


13.64 


45.18 


3.58 


49.96 


13.65 


50.04 


3.91 


55.07 


13.37 


54.93 


4.23 


60.17 


12.78 


59.83 


4.48 


65.26 


11.92 


64.74 


4.68 


70.30 


10.81 


69.70 


4.77 


75.32 


9.53 


74.68 


4.67 


80.29 


8.04 


79.71 


4.39 


85.23 


6.40 


84.77 


3.94 


90.17 


4.56 


89.85 


3.10 


95.10 


2.57 


94.90 


1.89 


100 


0 


100 


0 
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TABLE II.- CMP' a FOR I30LATBD AIRFOIL AND FOR AIRFOILS IH CASCADE HAllONAL ADWSORT 



COMMITTEE FOR AFRONAUTinS 



9 

(radlana! 


Isolated 


Cascade 


Ax 


Ay 


dAx 

dcp 


dAy 

dtp 


z 


k 


Ax 


A£ 
ir 


dAx 


dAy 


X 

IT 


k 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 


-0.125 

-.152 
-.156 
-.139 
-.105 
-.0515 
.0105 
.0744 
.122 
.146 
.148 
.127 
.0887 
.0350 
-.0121 
-.02Sfi 

-.0042 
.0223 
.0319 
.0317 
,0277 
.0063 
-.0294 
-.0788 


0.027 

.076 
.133 
.188 
.235 
.266 
.273 
.258 
.215 
.160 
.101 
.051 
.014 
.001 
.026 
067 
.092 
.092 
.074 
.057 
.032 
.009 
-.007 
-.006 


-0.137 

-.0655 
.0130 
.0985 
.175 
.228 
.244 
.228 
.150 
.0553 
-.0401 
-.116 
-.180 
-.215 
-.122 

.105 
.0693 
.005 
-.003 
-.0458 
-.105 
-.168 
-.206 


0.165 

.208 
.220 
.205 
.160 
.0770 
-.0250 
-.124 
-.196 
-.214 
-.200 
-.162 
-.107 
.0245 
.149 

.0497 
-.0596 
-.0682 
-.0646 
-.101 
-.0850 
-.0360 

.0721 


0.987 

.923 
.810 
.653 
.462 
.252 
.0314 
-.187 
-.403 
-.604 
-.776 
-.906 
-.982 
-.998 
-.936 
-.774 
-.529 
-.239 
.0528 
.335 
.594 
.799 
.937 
.996 


0.197 

.372 
.528 
.645 
.721 
.760 
.778 
.766 
.750 
.685 
.567 
.394 
.192 
.0654 
.412 
74S 
.963 
1.031 
1.005 
.964 
.829 
.616 
.347 
.0966 


-0.729 

-.610 
-.346 
-.190 
-.0856 
-.0123 
.0407 
.0776 
.113 
.122 
.146 
.0739 
-.0042 
.227 
.382 
.371 
.355 
.303 
.257 
.177 
.0866 
-.0484 
-.216 
-.490 1 


-0.155 

-.029 
.008 
.005 
-.003 
-.013 
-.026 
-.036 
-.049 
-.060 
-.074 
-.095 
-.047 
-.012 
-.071 
-.113 
-.149 
-.176 
-.202 
-.225 
-.247 
-.264 
-.273 
— .258 


-0.315 

1.028 
.800 
.466 
.320 
.230 
.170 
.122 
.084 
.0591 
.0281 
-.557 
.293 
1.087 

.10 
-.04 
-.12 
-.18 
-.235 
-.28 
-.395 
-.569 
-.794 
-1.218 


1.695 

1.024 
.0408 

-.0627 

-.110 

-.14 

-.14 

-.14 

-.14 

-.16 

-.28 
.0480 
.896 

-.304 

-.745 

-.46 

-.35 

-.30 

-.28 

-.27 

-.24 

-.14 
.03 
,631 


0.755 

.982 
.989 
.927 
.848 
.762 
.667 
.560 
.436 
.273 
.0530 
-.325 
-.864 
-.982 
-.853 
-.744 
-.638 
-.544 
-.448 
-.354 
-.247 
-.124 
.0499 
.316 


0.366 

.181 
.158 
.510 
.825 
1.095 
1^289 
l;392 
1.386 
1.264 
1.045 
.810 
.354 
.174 
-.503 
.796 
1.023 
1.175 
1.247 
1.256 
1.144 
.974 
.798 
.577 
















Interpolated 


0.5 
10.5 
11.5 
22.5 
23.5 














-.720 
.130 
.010 
-.340 
-.635 


-.085 
-.087 
-.082 
-.270 
-.219 


.65 
-.24 
-.36 
-1.01 
-1.16 


1.55 
-.27 
.62 
.15 
1.30 


.914 
-.107 
-.598 
.166 
.524 


.265 
.953 
.590 
.684 
.470 



TABLE III.- TRANSFORMATION CONSTANTS FOR TABLE II 



o 

Q 



H!Z! 
W 

□ < 
W H 
O CO 

ss o 





0 


(deg) 


K 


CPN 


Cprp 


Pt 


dcj 
da 


T 


R 


Cascade 


1 


45° 


0.258 


22°57' 


190°0» 


10°4' 


4.65 


0.197 








Isolated 


0 


0 


CO 


-10° 10' 


191°22« 


11°22' 


6.86 


.0209 


1.0912 



TABLE IV.- VARIATION OF PLOW ANGLES V/ITII ELADE LOADING 





a 


\ 


Ki 


M 




6n 


6ip 


Pu-Pd 


Cascade 


0 


-10°4» 


34°56' 


34°56' 


34°5G' 


0 


32°34' 


17°3« 


0 


.5 


-3°54» 


41°6' 


46°6' 


35° 19' 


10°47» 


21°24' 


17°26 » 


.329 


1.0 


2021' 


47°21« 


55°27' 


35°43' 


19°44' 


12°3' 


17°50' 


.734 


1.17 


4°31» 


49°31» 


68°23' 


35°51' 


22°32» 


9°7' 


17°58' 


.891 


(Ideal) 


















1.3 


6°10' 


51°10» 


60°24' 


35°54' 


24°30' 


7°6' 


18°1» 


1.013 


1.5 


8°45» 


53°45' 


63°27' 


36°12' 


27°15' 


4°3' 


18°19' 


1.214 


2.0 


15°25» 


60°25' 


70°11' 


36°52' 


33°19' 


-2°41' 


18°59 ' 


1.741 


Isolated 


0 


-11°22» 










33°52 ' 


15°45' 




.5 


-7°11» 










29°41' 


19°56 » 




1.0 


-2°59» 










25°29' 


24°8' 




1.28 


-0°36» 










23°6» 


26°31' 




(Ideal) 


















1.5 


1°16» 










21°14' 


28°23.» 




2.0 


5°36» 










16°54' 


32°45' 





!>-■ 



!2i 
O 



01 



to 



TABLE v.- CASCADE CIVP BY INVERSE METHOD 
First approximation 



a 

>^ 

M < 
to H 
O 

!2!0 







II 






Xl 
Tr 


^1 


{ci = 1.26) 


1 


0 


15 


-0.717 


-0.241 


-0.411 


1.673 


0.768 


0.344 


1.598 


1 Upper 


1 


_ cpT 


— . X±rt 


• fO 1 


1 npR 


QRQ 




1 P4fl 


J surface 


2 




_ "StWW 
~ •«3 1 f 


_ r>7Qfl 


7m 








7.'S7 

. f O r 




3 




■ COO 


_ nft^R 

— • wOwO 


Tip A 


— ■ wtst: w 


904 


599 


■ 803 




4 






— . UD 1 «J 






. OJ.VJ 




R1 Q 




5 










- 0.'^22 


712 

■ r Xm 


1 -194 

X • X 


829 




6 




_ -134 


- - 0902 


.0862 


.0342 


.607 


1.366 


.832 


^ Lower 
surface 


7 








.0694 


.0201 


.494 


1.430 


.034 

■ WW^ 


8 




_ 0073 


--0R16 


.0666 


.129 


.364 


1.359 

X ■ WW w 


.832 

• Www 


9 






-.0718 


.0933 


.0547 

■ \« w^ f 


.203 


1.228 

X • *•»*-• W 


.839 

• wWv 




10 






-.0566 


.1 86 


.248 


- .0158 


.979 

« V f w 


.840 




11 






-.0312 


213 


.176 


- .356 


.691 

• W vX 


.793 

• f i7W 




12 




.0747 


-.0125 

4 V r ill fir 1 / 


.575 


.502 


-.867 

• Ww f 


.319 

• W^v 


.708 

• f WW 




13 




.336 


.0021 


1 21 6 


- .418 


- 974 


.197 

■ Xv f 


. 614 




14 




.534 


••.0695 


.265 


_ .958 


- .826 

• WM W 


.587 

■ WW ■ 


.809 

• WW w 




15 




.545 


-.130 


.0172 

• W ^ ff M 


- .683 

■ www 


-.706 

• • wW 


.879 

■ W I V 


.916 

■ W JUW 




16 




.532 


-.183 


-.128 


-.533 


-.597 


1.073 


1.018 


^ Upper 
surface 


17 




.468 


-.223 


-.210 


-.441 


-.507 


1.191 


1.104 


18 




.409 


-.261 


-.287 


-.410 


-.414 


1.232 


1.184 


19 




.306 


-.294 


-.366 


-.373 


-.327 


1.194 


1.268 




20 




.188 


-.322 


-.464 


-.272 


-.230 


1.055 


1.386 




21 




.0358 


-.343 


-.624 


-.225 


-.112 


.938 


1.442 




22 




-.162 


-.357 


-.919 


.0363 


.0527 


.748 


1.533 




23 




-.460 


-.341 


-1.245 


.624 


.317 


.563 


1.551 , 





O 



TABLE v.- CASCADE CMP BY INVERSE JSETHOD - Concluded 

Second approximation 







TT 


rt At-— 

dtp 


dqp 


^2 

TT 


^2 


iCj - l.fio; 




0 


X X 
12 


-0.676 


-0.212 


-0.220 


1.269 


0.802 


0.334 


1.588 


L Upper 


1 


— . OriO 


_ T pf\ 
— • J.bW 


a ocw 




1 0 

X . W 


124 


1.508 


Isurfaoe 


2 




/ OA 




• ^OO 


OXK 


QR7 




77R 




3 




- .clU 




XDA 




. 090 


t^7(^ 


R41 
. 04X 




4 




- .£00 


r\A 1 


OCR 


m xn 




a yvJO 


R4P 




5 






- . U4Uy 


* J.O1C 


m no 
— « UXoo 


• 1 X / 


1 1 flfi 
X . xoo 






6 




- . ley 


r\A on 






.0x0 


1 xxo 


R4P 


Lower 


7 




- . UbUi 


- .Uoyo 


• uyco 




. i^XX 


1 XQ4. 


R44 


^aurfabe 

J 


8 




— « U&Uc 


- • UOoO 


nr? AT 
»u f o± 


nan n 


. ooo 


1 Xf.Q 


RXP 


9 






- .uom 


1 OA 




0X1 
. /SOX 


1 1 QR 

X . Xi70 


R4S 




10 






AO XX 


OT Q 


. XC U J 


01 Q7 
. UXO / 




R4R 




11 






— . yjGUO 


• X 1 o 




_ XI Q 
— • 0x9 


. wVJX 


-790 


> 


12 






— . \J\JO 1 I 


XRA 


xon 




.XI fl 

. 0x0 


RQ1 




13 




• OOO 






_ Atrn 
m ttJa 


_ 074. 


.188 
. xoo 


-657 


■ 


1 A 
±'± 




. 004 


— . W f 


nR-xv 

• v/Ow r 




— . OmO 


544 


862 




15 




■ OOc 


_ TZT 
— . 






- 71 fi 

— . 1 X\J 




996 

• V <7 V 




16 




.503 


-.179 


-.240 


-.408 


-.611 


.973 


1.104 




17 




.423 


-.214 


-.294 


-.376 


-.524 


1.092 


1.184 


Upper 
'^surface 


18 




.364 


-.247 


-.359 


-.345 


-.431 


1.135 


1.262 


19 




.250 


-.274 


-.418 


-.306 


-.348 


1.120 


1.327 


20 




.150 


-.296 


-.454 


-.241 


-.244 


1.079 


1.343 




21 




-.0010 


-.316 


-.660 


-.218 


-.126 


.908 


1.456 




22 




-.180 


-.329 


-.944 


.0611 


.0440 


.730 


1.531 




23 




-.488 


-.310 


-1.246 


.653 


.308 


.555 


1.528 


✓ 



TABLE VI.- ASSOCIATED CONSTANTS FOR CMP OP TABLE V 

1^ = 1, P = 45°] 



Approximation 


0 

(a) 


p 

(a) 


K 


Ik 


CI>ji 


Pt 

(a) 




T 


1 


1.001 


42°30' 


0.247 


21°10' 


189°46' 


7°5' 


4.69 


0.167 


2 


1.001 


41O3O' 


.237 


14°17» 


189°40» 


5°17« 


4.75 


.172 



Measured from "longeat-llne" chord of airfoil. 

TABLE VII.- CASCADE PLOW ANGLES FOR DERIVED AIRFOIL SECTIONS 



Approximation 


H 


a 

(a) 


\ 


\i 


^6. 


A\ 


^N 




-Pu" ^d 


1 


1.26 


8^30 » 


51° 


60°4' 


36°18» 


23°46 » 


6°32' 


-7°53' 


0.98 


2 


1.26 


10°6' 


51°36' 


60°31' 


37°4» 


23°27' 


-8°42' 


-7°56' 


.99 



Measured from "longest-line" chord of airfoil. 
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Figure ^— Cascade f /at -plate transformaf/on. 
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Figure J. — Source - y^ortex system 

for cascade f/Okv. 



Given airfoil /n cascade 



V Zero approximation ^ upper surface 
A Z^rc a/:?/proy//^a/'/o/7 ^ lower s6/rf^7^:ff 
-Or -Firs f £?/:?/?rcvr//??a//j>^ NATIONAL ADVISORY 

□ 5econd d7/?yOrCir//rxp//a/7 COMMITTEE FOR AEBONAUTJCS 

A Third £7yO/prajc//77(2//a/? _ 
d Interpolated points 




Loiver surface 



Figure 4. — Direct CAfF mefttod for a /r foil In cascade. s=4J° , 



Fig. 6 
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Upper surface 



A .G 3 . /.O 

Chord 

r/gure e— l^e/oc/fy d/sfr/£)uf/ons af- c^-O. 
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Fig. 7 




Figure 7. — ^^e/oc/fij d/sfr/butions ^ /50/afed airfoil, 




F/gureS. — l/^/oc/fy d/'sfr/'/bu/'/on^ ^a/rfo// m casccfcfe^ 
^^^^^^ 0'=/, 
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Fig. 9 



— Prescribed i/eJocjfy 

- - Inifia/ we/oc/ty^ c^=^/.3 
o First approximation 
B Second ^yo/?r€>x/yn 4?7Vo/? 



Upper surface 




lower surface 



Figure 9— /eiocify distriJbutions in cascade 

Jby inverse method. 



Ifiifta/ airfoif 




Hgure 10. — -Cascade airfoih cfer/ved Jt^y /nt^er-se method. 
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Fig. 11 




Figure //. — CMF denVaftVes by inrerse method for 
a/rfoifs in cascade. 



